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Abstract 

><" 

; QCD is constructed as a lattice gauge theory in which the elements 

of the link matrices are represented by non-commuting operators acting 
in a Hilbert space. The resulting quantum link model for QCD is for- 
mulated with a fifth Euclidean dimension, whose extent resembles the 
inverse gauge coupling of the resulting four-dimensional theory after di- 
mensional reduction. The inclusion of quarks is natural in Shamir's vari- 
ant of Kaplan's fermion method, which does not require fine-tuning to 
approach the chiral limit. A rishon representation in terms of fermionic 
constituents of the gluons is derived and the quantum link Hamiltonian 
for QCD with a U{N) gauge symmetry is expressed in terms of glueball, 
meson and constituent quark operators. The new formulation of QCD 
is promising both from an analytic and from a computational point of 
view. 
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1 Introduction 



Solving QCD is among the most challenging problems in theoretical physics. A 
non-perturbative formulation of QCD is provided by Wilson's lattice gauge theory 
which maps the problem to one of classical statistical mechanics. Over the past 
twenty years a variety of tools have been developed to solve lattice field theories. At 
present the most powerful tool is the Monte Carlo simulation of the partition func- 
tion of the corresponding classical statistical mechanics system. However, the most 
efficient numerical algorithms for solving lattice QCD suffer from critical slowing 
down when the continuum limit is approached and thus exhaust even the biggest 
supercomputers. The use of improved or even "perfect" actions may alleviate this 
problem 0, because it may allow one to extract continuum physics from rather 
coarse lattices. On the other hand, to obtain high-precision numerical data us- 
ing the standard numerical techniques will still require large resources of computer 
power. Hence, it is certainly reasonable to look for new formulations of the QCD 
problem. 

Here we propose an alternative non-perturbative approach to QCD in the frame- 
work of quantum link models. Such models were first discussed by Horn 0, and 
studied in more detail by Orland and Rohrlich |^. In these models the classical sta- 
tistical mechanics problem of standard lattice gauge theory is replaced by a problem 
of quantum statistical mechanics. In particular, the classical Euclidean action is re- 
placed by a Hamilton operator. As a consequence, the elements of the link matrices 
that are ordinary c-numbers in the standard formulation of lattice gauge theory, 
turn into non-commuting operators acting in a Hilbert space. Recently, quantum 
link models have been formulated with a fifth Euclidean direction and it has been 
shown how they are related to ordinary 4-d Yang-Mills theories via dimensional 
reduction [Q. The previous work on quantum link models was limited to U{N) 
and SU{2) gauge groups and it was not clear how to generalize the construction 
to SU{N). The quantum link models constructed in this paper have an SU{N) 
gauge symmetry and provide a new non-perturbative formulation of QCD. We hope 
that this formulation will lead to alternative ways for attacking this long-standing 
problem, both with analytic and with numerical methods. 

Quantum link QCD is quite different from the standard formulation of lattice 
QCD. The main novel feature of quantum link models is that the gluonic Hilbert 
space on each link is finite. While the Hilbert space in the Hamiltonian formula- 
tion of standard lattice QCD contains all representations of SU{3) on a link, the 
corresponding Hilbert space of an SU{N) quantum link model consists of a single 
representation of SU{2N). This is achieved by formulating the theory in 4 -|- 1 di- 
mensions. In order to represent full QCD, of course, quarks must also be included. 
Even in the standard formulation of lattice gauge theory, the Hilbert space of the 
quarks is finite and indeed the representation of quarks in quantum link QCD is not 
very different. However, the boundary conditions of the quarks in the fifth Euclidean 
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direction must be carefully adjusted in order to ensure their proper dimensional re- 
duction. Antiperiodic boundary conditions do not work, because they do not lead to 
a long enough fermionic correlation length in four dimensions. Periodic boundary 
conditions, on the other hand, allow an arbitrarily long 4-d fermionic correlation 
length. The resulting 4-d fermion resembles in many respects the standard lattice 
Wilson fermion. In particular, with periodic boundary conditions for the quarks in 
the fifth direction, quantum link QCD suffers from a fine-tuning problem, when one 
wants to approach the chiral limit. Fortunately, we can do better than that. 

Kaplan proposed to represent 4-d chiral fermions by working with a domain 
wall in a five- dimensional world 0. In his construction the fifth dimension arises 
for completely different reasons than in quantum link models. However, since we 
already have five dimensions, it is natural to use them also to protect the chiral 
symmetries of the fermions. The scheme that fits most naturally with quantum 
link models is Shamir's variant [0] of Kaplan's proposal, which works with a 5- 
d slab (which is finite in the fifth direction) with open boundary conditions for 
the fermions at the two boundaries. Due to the boundary conditions, a left-handed 
fermion is bound to one side of the slab, while a right-handed fermion appears on the 
opposite side. Remarkably, the correlation length of these fermions is exponentially 
large in the depth of the slab. This is exactly what one needs for them to survive 
the dimensional reduction in quantum link models. The corresponding geometry 
is shown in fig.l. When Shamir's boundary conditions are used for the quarks in 
quantum link QCD, no fine-tuning is necessary to approach the chiral limit of the 
theory. This is essential, because recovering the chiral properties of the quarks is one 
of the main problems in numerical simulations of Wilson fermions. Recently, using 
the standard formulation of lattice gauge theory, it was demonstrated convincingly 
that Shamir's proposal is indeed practical from a computational point of view 
and that it leads to a great improvement of the chiral properties of quarks on the 
lattice. 

It is reassuring that the formulation of quantum link models in five dimensions 
fits naturally with this elegant solution of the fermion doubling problem in vector-like 
theories including QCD. Should similar ideas ultimately lead to a lattice construction 
of chiral gauge theories, one may expect that the theory can also be formulated as a 
quantum link model. Once the lattice chiral fermion problem is solved, the door may 
be open to a non-perturbative lattice formulation of supersymmetry. In that case 
quantum link models may offer the most natural framework, because they treat the 
Hilbert spaces of bosons and fermions on an equal footing. In particular, in quantum 
link models the bosonic Hilbert space is also finite (on a finite lattice). 

Due to the fact that the gluonic Hilbert space of quantum link QCD is finite, one 
can reformulate the theory in terms of fermionic constituents of the gluons, which 
we call rishon^ following Jewish tradition |p. The rishons carry color and transform 

^Rishon is Hebrew and means "first". 
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Figure 1: Dimensional reduction in slab geometry: The gluonic correlation length 
1/m is exponential in j3. Left and right-handed quarks live on the two sides of 
the slab. The overlap of their wave functions and induces a mass fi. The 
corresponding correlation length is exponentially larger than 1/m. 



in the fundamental representation of SU{N). Each quantum link variable can be 
expressed as a rishon-anti-rishon pair. The rishon dynamics is particularly simple. 
The Hamilton operator describes the hopping of rishons from one end of a link to the 
other, as illustrated in fig.2. The rishon formulation of quantum link QCD resembles 
the Schwinger boson and constraint fermion constructions of quantum spin systems 
[|T0[| and provides theoretical insight, that may lead to new analytic approaches to 
QCD, for example, in the large N limit. The Hamilton operator of quantum link 
QCD with a U{N) gauge group can also be expressed in terms of glueball, meson and 
constituent quark operators. These objects consist of two rishons, two quarks and a 
quark-rishon pair, respectively. It is conceivable that the theory can be bosonized in 
terms of these degrees of freedom. If large techniques can be applied successfully, 
major progress in the non-perturbative solution of QCD and other interesting field 
theories is to be expected. It is intriguing that U{N) quantum hnk QCD has two 
representations — one exclusively in terms of color-triplet quark and rishon fermions. 
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Figure 2: Rishon dynamics: The trace part of the Hamiltonian induces a hopping 
of rishons of various colors around a plaquette. The determinant part shifts a color- 
neutral rishon baryon from one end of a link to the other. 



the other in terms of color-singlet bosons. 



In contrast to standard lattice gauge theories, the configuration space of quan- 
tum link models is discrete and their path integral representation resembles that of 
quantum spin models. Hence, it seems plausible that highly efficient cluster algo- 
rithms similar to those for quantum spins become available for quantum link 
models. For classical spin models the development of cluster algorithms [|12| , |13[ gave 
rise to a tremendous improvement of numerical data and has led to a satisfactory 
high-precision numerical solution of various models. Despite many attempts, it has 
so far not been possible to construct efficient cluster algorithms in the standard for- 
mulation of lattice QCD. The 2-d classical 0(3) model is in many respects similar 
to 4-d Yang-Mills theories. The most efficient way to simulate the 0(3) model is 



provided by the Wolff cluster algorithm 



On the other hand, the same model 



can also be obtained from the 2-d antiferromagnetic quantum Heisenberg model via 
dimensional reduction |T1|. For the Heisenberg model a very efficient loop cluster 
algorithm is available [|1^, |T^, which has recently been modified to work directly 
in the Euclidean time continuum 
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This algorithm has been used to simulate 
the correlation length of the 2-d classical 0(3) model that results from dimensional 
reduction [|l^ . The quality of the numerical data is compatible with those obtained 



from the Wolff algorithm applied directly to the 2-d classical model. For example. 
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using a finite size scaling technique, correlation lengths of up to about 37000 lattice 
spacings could be handled successfully. This suggests that, if a cluster algorithm 
can be constructed for quantum link QCD, its numerical simulation should be more 
efficient than that of standard lattice QCD, for which no cluster algorithm has been 
found. 

At present, it is unclear if superior numerical methods can be developed in the 
quantum link formulation of lattice QCD. However, this seems possible, because, for 
example, for a U{1) quantum link model a cluster algorithm operating in the con- 
tinuum of the fifth Euclidean direction has already been constructed [|T^ . Quantum 



link models with quarks resemble standard models of quantum statistical mechanics, 
as, for example, the Hubbard model. This may eventually lead to problems with 
their numerical simulation, because such models are known to suffer from the noto- 
rious fermion sign problem. However, based on a fermion cluster algorithm [|19[, a 



recently proposed method may circumvent this difficulty. However, we postpone 
further considerations of numerical simulations of quantum link QCD with dynami- 
cal quarks. The main issue of this paper is to present the quantum link formulation 
of QCD. 

The paper is organized as follows. In section 2 quantum link models with U{N) 
and SU{N) gauge symmetries are constructed and in section 3 it is discussed how 
these models are related to ordinary 4-d Yang-Mills theories via dimensional reduc- 
tion. In section 4 quantum link models are reformulated in terms of rishons — the 
fermionic constituents of the gluons. The inclusion of quarks is discussed in sec- 
tion 5 within the framework of Shamir's variant of Kaplan's fermion method. The 
representation of the U{N) gauge invariant quantum link Hamiltonian, in terms of 
glueball, meson and constituent quark operators, is worked out in section 6. Finally, 
section 7 contains our conclusions. 



2 U{N) and SU{N) Quantum Link Models 

Quantum link models were first discussed by Horn P]. He constructed models with 
U{N) and SU{2) gauge groups, but he pointed out that his construction can not 
be generalized to SU{N). Further, the recent application of quantum link models 
to Yang-Mills theory via dimensional reduction was discussed only for SU{2) 
Here we explicitly construct quantum link models with an SU{N) gauge symmetry. 
Obviously, this is essential for the quantum link formulation of QCD with an SU{3) 
gauge group. 

Let us first recall the standard formulation of lattice gauge theory. In that case 
there is an SU{N) matrix u^^^ associated with each link {x,fi) on a 4-d hypercubic 
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lattice. The action is given by 

S[u] ^- J2 TrK^iix+A,i^'"l+i>,^4,J> (2-1) 

where the dagger denotes Hermitean conjugation. The action is invariant under 
SU {N) gauge transformations 

u'x,n = exp(ia^ • X)u^^^ exp(-ia^+^ • A), (2.2) 

where A is the vector of Hermitean generators of SU{N) with the usual commutation 
relations 

[A«, A^] = 2i/„,eA^ TrA"A^ = 25"^ (2.3) 
The path integral is given by 

[ Vu&cp{-\s[u]), (2.4) 

where g is the non-Abclian gauge coupling. Formally, the above system can be 
viewed as one of classical statistical mechanics. Then the action plays the role of 
the classical Hamilton function and g'^ plays the role of the temperature. The 4-d 
SU (N) lattice gauge theory is believed to have only one phase, in which the gluons 
are confined. Due to asymptotic freedom, the continuum limit of the lattice model 
corresponds to gi ^ 0. 

Let us now construct a quantum link representation of SU{N) lattice gauge 
theory. We replace the classical action (or Hamilton function in the language of 
statistical mechanics) by a quantum Hamilton operator 

Here the elements of the N x N link matrices Ux,^ are operators acting in a Hilbert 
space — not just c-numbers like in the standard formulation. Naturally, the dagger 
now represents Hermitean conjugation in both the Hilbert space and the N x N 
matrix space. Note, however, that the trace in the above expression is only over 
the N X N matrix space and not over the Hilbert space. Gauge invariance of the 
quantum link model requires that the above Hamilton operator commutes with the 
generators Gx of infinitesimal gauge transformations at each lattice site x, which 
obey the standard local algebra 

[G:,Gl]^2i5,yUcGl. (2.6) 

Gauge covariance of a quantum link variable follows by construction if 

^ Ylewi-i^y ■ Gy)Ua:,,j,'[lex.p{ia:, ■ G^) = exp(ia^ • A)^/^,,^ exp(-ia3,+^ • A), 
y z 

(2.7) 
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where Ox Gxp{idx ■ Gx) is the unitary operator that represents a general gauge trans- 
formation in Hilbert space. The above equation implies the following commutation 
relations 

[Gx, Uy^^] = Sx^y+fiUy^^X — Sx^yXUy^f^. (2.8) 

In order to satisfy these relations we introduce 

Gx = "^{Rx-fi,!! + Lx,^), (2.9) 

where Rx^fi and Lx^fj. are generators of right and left gauge transformations of the 
link variable Ux,^- Suppressing the link index (x,/i), their commutation relations 
take the form 

R'] = 2tUcR', L'] = 2tU,L% L'] = 0, (2.10) 

i.e. R and L generate an SU{N)fi^® SU{N)l algebra on each link. The R and L op- 
erators associated with different links commute with each other. The commutation 
relations of eg. (p78|) imply 



[R, U] = UX, [L, U] = -\U. (2.11) 

The above commutation relations are identical to those of the standard Hamiltonian 
formulation of lattice gauge theories. There they are the canonical commutation re- 
lations between link variables, which play the role of coordinates, and electric field 
variables R and L, which play the role of canonical conjugate momenta. In the stan- 
dard formulation the commutation relations are realized in an infinite dimensional 
Hilbert space of square integrable wave functionals of the link variables. This Hilbert 
space is necessarily infinite dimensional, because one insists that the elements of a 
link matrix commute with each other. Quantum link models, on the other hand, 
have a finite dimensional Hilbert space for each link, by allowing the elements of the 
link fields to have non-zero commutators. It is important to realize that this does 
not destroy the gauge symmetry of the problem. 

For each link the above relations can be realized by using the generators of 
an SU{2N) algebra, with the SU{N)l ® SU{N)r algebra embedded in it. In the 
fundamental representation of SU{2N), for example, we can write 

R = 





Uij = Rellij + i ImUij, {U'^)ij = Rellji — i ImU. 



The M^*-'^ are a set of N'^ matrices (one for each index pair (ij)) of size N x N with 
^kP ~ ^ii^jk- The commutators of real and imaginary parts of the elements of the 
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link matrices take the form 

[ReUij, ReUki] = [Imf/jj, ImUki] = -i{5ik ImA^-i ■ R + 5ji ImA^fc ■ L), 

^ ^ ^ ^ 2 

[Ref/ij, Imf/fci] = i{5ik ReXji ■ R - Sji ReXik ■ L + —6ikSjiT). (2.13) 

In particular, [Uij,Uki] = [{U'^)ij, {U'')ki] = 0. Again, we should emphasize that 
the commutation relations are local, namely all commutators between operators 
assigned to different links are zero. 

The additional generator T is given by 

The real and imaginary parts of the N"^ matrix elements Uij are represented by 2N'^ 
Hermitean generators of SU{2N). Together with the 2{N'^ — 1) generators R and L 
of right and left gauge transformations and the generator T, these are the 

2N'^ + 2{N^ - 1) + 1 = AN"^ - 1 (2.15) 

generators of SU{2N). The generator T commutes with the generators of SU{N) 
gauge transformations, i.e. 

[T,R] = [T,L] = 0, (2.16) 
but not with the elements of the link matrices, because 

[T, U] = 2U, [T, U^] = -2Ul (2.17) 

This relation implies that 

= 9 5I(^^-A,M ~ ^a;,/i) (2.18) 

generates an additional U{1) gauge transformation, i.e. 

= Y[exp{-iayGy)Ux,^,Ylexp{ia^G^) = exp(ia^)f/^^^ exp(-m^+^). (2.19) 
y z 

Indeed the Hamilton operator of eq. (p.5|) is also invariant under the extra U{1) 
gauge transformations and thus describes a U (N) lattice gauge theory. The SU (2) 
quantum link model of ref.[^ works with the generators of SO (5) instead of SU{A). 
That construction does not have the additional U{1) gauge symmetry, but it can 
not be generalized to SU {N) as will be explained below. 

The question then arises how quantum link models can be used to represent QCD. 
The answer is surprisingly simple. All one needs to do is to break the additional 
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f/(l) gauge symmetry by adding the real part of the determinant of each hnk matrix 
to the Hamihon operator, such that now 

H = J Tr[f/.,^t/.+^,,f/i+,,^t/t_J + J'Y: [detf/.,^ + detf/J^^]. (2.20) 

Since all elements of Ux,^ commute with each other (although the real and imagi- 
nary parts of an individual element do not commute) the definition of detf/j,^ does 
not suffer from operator ordering ambiguities. By construction the above Hamilto- 
nian is invariant under SU{N) but not under the extra U{1) gauge transformations. 
However, there is a subtlety that needs to be discussed. In the fundamental repre- 
sentation of SU{2N), that was discussed above, the operator that represents detUx^fj, 
turns out to be zero. Hence, in that case even the Hamiltonian of eq. (|2.20|) has a 
U{N) gauge invariance. On the other hand, since the above commutation relations 
can be realized with any representation of SU{2N), we can use higher represen- 
tations in which the determinant in general does not vanish. In section 4 it will 
become clear that the (2iV)!/(iV!) ^-dimensional representation of SU{2N) is the 
smallest with a non- vanishing determinant. For QCD the Hilbert space of the cor- 
responding SU{3) quantum link model then is the direct product of 20- dimensional 
Hilbert spaces associated with each link. Specifying the state on a link thus requires 
just 5 bits (because 2"^ < 20 < 2^). In the standard formulation of lattice gauge 
theory, on the other hand, in the computer an SU{3) link variable is usually rep- 
resented by 18 real numbers, which, for example, on a CRAY corresponds to 912 
bits. 

At this point it is easy to see why there is a special construction for SU (2) us- 
ing a 4-dimensional representation |^. Under gauge transformations the quantum 
link operator U transforms as an {A^, A^} tensor. In SU (2) the fundamental rep- 
resentation {2} and its conjugate {2} are unitarily equivalent. Consequently, one 
can construct another quantum link operator using f/^ with the same SU{2) gauge 
transformation properties as U . A linear combination of the two operators leads to 
a quantum Hamiltonian with SU{2) — but not U{2) — gauge symmetry. Indeed, 
adding them with equal weight gives the 5*0(5) construction of 

It should be clear that the Hamiltonians of eqs.(|2.5|) and ( f^.20D represent just 
the simplest members of a large class of U{N) and SU{N) invariant quantum link 
models. Like in ordinary lattice gauge theories, the action (in this case the Hamilto- 
nian) can be improved by including more complicated terms, for example, a six-link 
loop around a double plaquette. Further, it is possible to break the extra U{1) 
gauge symmetry by antisymmetrized combinations of various paths connecting two 
lattice points — not just by the determinant operator on a single link. In that 
way one can construct SU{N) invariant quantum link Hamiltonians even with the 
fundamental representation of SU{2N), but with more complicated interactions. In 
the 15-dimensional representation of SU{6) one can construct a more complicated 
plaquette action that leads to an SU (3) invariant theory, even without the link de- 
terminant term. All these options may become important when the question of the 
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most efficient numerical treatment of quantum link models will be addressed. 

For the moment we limit ourselves to the discussion of one natural extension of 
the Hamilton operators from before. As it was discussed in ref. 0, the operators 
Tx,ii-, Rx,fi and Lx^fj, resemble Abelian and non-Abelian 5-d electric fluxes associated 
with the link (x, fi). On the other hand, from a 5-d point of view the plaquette term 
in the above Hamilton operators represents the magnetic part of the energy only. 
Hence, it would be natural to add a 5-d electric term as well. This corresponds to 

H = J ^ Tr[f/.,^f/.+^,,t/i+,XJ + '^T [dett/x,;. + det[/i^^ 

+ J"T.[Rl, + LlJ + r'T. Tl- (2.21) 

X,fJ, X,fJ, 

This Hamiltonian is still gauge invariant, because -R^^^, Ll^^ and T^^ are Casimir 
operators of SU{N)ji ® SU{N)l ® f/(l), which commute with the generators of 
all gauge transformations. In the fundamental representation of SU{2N) the above 
combination of Casimir operators is proportional to the unit matrix and has no effect 
on the dynamics. This is not the case for the higher-dimensional representations of 
SU{2N). In principle, one can also include the other Casimir operators of SU{N)ji^ 
SU{N)l ® U{1) in the Hamiltonian. The above choice of Casimir operators is most 
natural, because it represents Abelian and non-Abelian 5-d electric field energies. In 
the following we present universality arguments, which suggest that quantum link 
models defined with a 4-d Hamiltonian (describing the evolution of the system in 
a fifth Euclidean direction) get dimensionally reduced to ordinary 4-d Yang-Mills 
theories. These arguments are based on symmetry considerations and they apply 
to all SU{N) invariant Hamiltonians discussed in this section. In the standard 
formulation of lattice gauge theories, one must include the electric terms in the 
action in order to get non-trivial dynamics. Since the real and imaginary parts 
of the elements of a quantum link operator do not commute, the dynamics of a 
quantum link model is already non-trivial without the electric terms. For simplicity, 
we restrict ourselves to J" = J'" = for the rest of this paper. 



3 Reduction from Five to Four Dimensions 

As discussed in detail in ref. quantum link models in 4-|- 1 dimensions are related 
to ordinary 4-d gauge theories via dimensional reduction. The Hamilton operator 
of the quantum link model, which is defined on a 4-d lattice, describes the evolu- 
tion of the system in a fifth Euclidean direction. The crucial observation is that 
a genuine 5-d non-Abelian gauge theory has deconfined massless gluons and thus 
an infinite correlation length. When periodic boundary conditions are imposed in 
the fifth direction and its extent /3 is made finite, the extent of the extra dimension 
is thus negligible compared to the correlation length. Hence, the theory appears 
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to be dimensionally reduced to four dimensions. Of course, in four dimensions the 
confinement hypothesis suggests that gluons are no longer massless. Indeed, as it 
was argued in ref . |^ , a glueball mass 

is expected to be generated non-perturbatively. Here e is the renormahzed dimen- 
sionful gauge couphng of the 5-d gauge theory. For large P the gauge coupling of 
the dimensionally reduced 4-d theory is given by 

1/g' = P/e\ (3.2) 

Thus the continuum limit (7 — of the 4-d theory is approached when one sends the 
extent f3 of the fifth direction to infinity. Hence, in contrast to naive expectations, 
dimensional reduction occurs when the extent of the fifth direction becomes large. 
This is due to asymptotic freedom, which implies that the correlation length ^ = 1 /m 
grows exponentially with jS. 

When a gauge theory is dimensionally reduced, usually the Polyakov loop in the 
extra dimension appears as an adjoint scalar field. Here we want to obtain pure 4-d 
Yang-Mills theory (without charged scalars) after dimensional reduction. This can 
be achieved if one does not impose Gauss' law for the states propagating in the fifth 
dimension, because the Polyakov loop is a Lagrange multiplier field that enforces 
the Gauss law. Formally, this can be realized simply by putting the fifth component 
of the gauge potential to zero, i.e. 

A, = 0. (3.3) 

In the infinite volume limit of the 5-d theory this restriction has no effect on the 
dynamics. With finite extent in the fifth direction, however, it deviates from the 
standard formulation of gauge theories. The leading terms in the low energy effective 
action of the 5-d gauge theory corresponding to the quantum link model take the 
form 

S[A,] = f^dx^ J d^x ^[Tr F^^F,, + ^Tr d^A^O.A^]. (3.4) 

We expect that the QCD quantum link model leads to a 5-d gauge theory charac- 
terized by the "velocity of light" c. Note that here /x runs over 4-d indices only. At 
finite (3 the above theory has only a 4-d gauge invariance, because we have fixed 
= 0, i.e. we have not imposed the Gauss law. On the other hand, for (3 = 00 
a full 5-d gauge symmetry is recovered, although the above action then still is in 
As = gauge. Since we are interested in dimensional reduction, a 4-d gauge sym- 
metry is sufficient for our purposes. On the level of the quantum link model, not 
imposing Gauss' law is achieved by simply writing the quantum statistical partition 
function as 

Z = Trexp(-/?i7). (3.5) 
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In contrast to the standard formulation of gauge theories we have not included a 
projection operator on gauge invariant states, i.e. gauge variant states also propagate 
in the fifth direction. The Hamilton operator of the quantum link model is defined 
on a 4-d space-time lattice and describes the evolution of the system in the fifth 
unphysical direction. In particular, all the information about the physical spectrum 
of the 4-d theory is contained in correlation functions in the Euclidean time direction, 
which is part of the 4-d lattice. Note that the physical Gauss law is properly imposed 
because the model does contain non-trivial Polyakov loops in the Euclidean time 
direction. In the continuum limit g ^ 0, which we approach by increasing the extent 
(3 of the fifth dimension, we are probing the low lying states in the spectrum of the 
4-d Hamilton operator of the quantum link model. The space-time correlations in 
these unphysical states of the 4-d Hamiltonian contain the information about the 
physical spectrum. 

It is useful to think of the dimensionally reduced 4-d theory as a lattice theory 
with lattice spacing /3c (which has nothing to do with the lattice spacing of the 
quantum link model). In fact, one can imagine performing a block renormalization 
group transformation that averages the 5-d field over cubic blocks of size f3 in the 
fifth direction and of size (3c in the four physical space-time directions. The block 
centers then form a 4-d space-time lattice of spacing (5c and the effective theory of 
the block averaged 5-d field is indeed a 4-d lattice theory. This is illustrated in fig. 3. 
A similar argument was first used by Hasenfratz and Niedermayer in their study 
of 2-d quantum antiferromagnets. 



The partition function of eq.( |3.5| ) can be written as a 5-d path integral of dis- 
crete variables — in the SU (3) case the states of the 20-dimensional representation 
of SU{6) on each link. In many respects this path integral resembles that of quantum 
spin systems, which can be simulated by very efficient loop cluster algorithms. Due 
to the discrete nature of the Hilbert space, one can even work directly in the contin- 
uum for the extra Euclidean direction [^. Indeed for a f/(l) quantum link model a 



cluster algorithm operating in the continuum has already been constructed [0 . It is 
plausible that cluster algorithms can also be constructed for non-Abelian quantum 
link models, which would allow high-precision simulations of gauge theories. 



4 The Rishon Representation of Quantum Link 
Models 

In this section we reformulate quantum link models using anticommuting operators 
describing fermionic constituents of the gluons. In contrast to composite models , 
our rishons are not physical particles, because they are absolutely tightly bound into 
gluons and there is no kinetic term for the rishons. Still, the rishons are convenient 
mathematical objects that provide theoretical insight into the physics of quantum 
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4-d ordinary lattice 
gauge theory 




Figure 3: Dimensional reduction to a 4-d lattice theory: A renormalization group 
transformation that averages the 5-d field over blocks of size (3 in the fifth direction 
and /9c in the four space-time directions results in an effective 4-d lattice gauge 
theory with lattice spacing (5c. 

link models. The rishon representation is interesting both from an analytic and from 
a computational point of view. First, it may lead to a new way to attack the large 
limit of QCD and other non-perturbative quantum field theories. In addition, 
the rishon formulation may also be useful in the construction of cluster algorithms. 

The algebraic structure of a quantum link model is determined by the com- 
mutation relations derived in section 2. The Hilbert space is a direct product of 
representations oi SU{2N) on each link, with the generators on different links com- 
muting with each other. We can thus limit ourselves to a single link, for which the 
commutation relations are 

R'] = 2tUcR', L'] = 2^/„fe,L^ L'] = [T, R-] = [T, L^] = 0, 
Uij] — UikXlj, [L°-, Uij] = —Xl^Ukj-, [T, Uij] = 2Uij, 
[ReUij,ReUki] = [ImUij,lmUki] = ImX^R'' + 5ji ImA^^L"), 

[Ret/,,-, ImUki] = ReA^^^i?" - 6ji ReX^^L^ + ^AkSjiT). (4.1) 
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These relations can be realized in a representation using anticommuting rishon op- 
erators c^ _|_^, c^|±^ with color index i G {1, 2, A^}. The rishon operators live at 
the left and right ends of the links and are characterized by a lattice point x and a 
link direction ±/i. They obey canonical anticommutation relations 

Under SU{N) gauge transformations the operators c and transform in the funda- 
mental and anti-fundamental representation, respectively. It is straightforward to 
show that the commutation relations of eq.( [4.1|) are satisfied when we write 

pa _ \^ „«t \a J J-a _ \^ Jt \a J 

id id 

i 

{Ux,^)ij — Cx,+fi'^x+fi,-iiy i^x,fi)ij ~ '^x+fi,-fi'^x,+fi- (4-3) 

Actually, the commutation relations are satisfied also when the rishons are quantized 
as bosons. It is crucial to realize that all operators introduced so far (including the 
Hamiltonian) commute with the rishon number operator 

■A^j/i = X^(cJ+^^_^C3.^^^_^ + Ca!+;iC^,+^) (4-4) 
i 

on each individual link. Hence, we can limit ourselves to superselection sectors of 
fixed rishon number for each link. This is equivalent to working in a given irreducible 
representation of SU{2N). Let us use the rishon representation to take a closer look 
at the determinant operator that we used to break the U {N) gauge symmetry down 
to SU{N). We have 



1 



detUx^jj, . £iii2...iff{Ux,^)iii'^ iUx,fi)i2i'2 ■■■{Ux,^)i 



ii ^i"^ i'2 ^2"^ ijv ^A^"^ //I c ^ 

Due to the antisymmetry of the e-tensors this operator would vanish for bosonic 
rishons. For fermionic rishons, on the other hand, we have 

^^^'-^x,fi ■ ^x,+fi^x+/i,-fi'^x,+fi'^x+fj.,-t^---'^x,+fi^x+fi,-fi- \^-^) 

Only when this operator acts on a state with exactly N" = N rishons (all of a dif- 
ferent color), it can give a non-zero contribution. In all other cases the determinant 
vanishes. This means that we can reduce the symmetry from U{N) to SU{N) via 
the determinant only when we work with exactly JV = N fermionic rishons on each 
link. The number of fermion states per link then is 

(2N\i2N)l 
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This is the dimension of the SU{2N) representation with a totally antisymmetric 
Young tableau with boxes (arranged in a single column). 



In the rishon representation the Hamilton operator takes the form 

x,fj,y^u i,j,k,m 



^ ^ • V-'x,+^i^x+fi,—^i^x,+ii^x+n,-fi---^x,+ii^x+n,-iJ, 

X,fJ, 

^X+fl,~fJ.^X,+fl^X+IJ.,~fl^X,+fl- ■ -^X+fJ^—fl^Xj+^l 

+ ^'11 [det c^,+/.4+^^_^ + det c^+/i,_^4,+^]. (4.8) 

The terms in brackets in eq. ( [4.8| ) represent color neutral glueballs formed by two ris- 
hons located at the same corner of a plaquette. The plaquette part of the Hamilton 
operator shifts single rishons from one end of a link to the other and simultane- 
ously changes their color. The determinant part, on the other hand, shifts an entire 
rishon-baryon — a color neutral combination of rishons — along the link. This is 
illustrated in fig.2. The rather simple rishon dynamics may facilitate new analytic 
approaches to QCD and may also be useful in numerical simulations of quantum 
link models. It is remarkable that the Hamilton operator can be expressed in terms 
of color neutral glueballs and rishon-baryons. This observation suggests a new ap- 
proach to the large limit, which is presently under investigation. 



5 Quantum Link Models with Quarks 

To represent full QCD, it is essential to formulate quantum link models with quarks. 
This is more or less straightforward, although some subtleties arise related to the 
dimensional reduction of fermions. Before we discuss the quantum link formulation 
of full QCD, let us review the standard formulation of lattice gauge theory with 
fermions. Here we concentrate on Wilson's method, although it is straightforward 
to construct quantum link models with staggered fermions along the same lines. 
The standard Wilson action with quarks is given by 

X,fJ,J^U 

+ [^xl^U^,p,iJx+f, - ^x+f^lt^ul^^lp:,] + M ^ ^p^i)^ 

+ 2 ^ [2^xV'x - ^xUx,i,iJx+f, - il)x+[Lul^f,'4'x]- (5.1) 
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Here ijjx and ipx are independent Grassmann valued spinors associated with the 
lattice site x, 7^ are Dirac matrices and M is the bare quark mass. The term pro- 
portional to r is the Wilson term that removes unwanted lattice fermion doublers at 
the expense of explicitly breaking chiral symmetry. In order to reach the continuum 
limit the bare mass M must be tuned appropriately. 

The guiding principle in the formulation of quantum link models is to replace the 
classical action of the standard formulation by a Hamilton operator that describes 
the evolution of the system in a fifth Euclidean direction. For the quarks we must 
also replace %jjx by ^^^75. Hence, the full QCD quantum link Hamiltonian is given 

by 

H = J ^ Tr[f/.,^f/.+^,,^i+,XJ + '^'E [detf/.,^ + dett/tj 

^ X,fl X 

+ ^ E [2^175^. - ^h5t/x,,^.+A - ^I+A75^i,^x]- (5.2) 

Here and '^x are quark creation and annihilation operators with canonical anti- 
commutation relations 

{VI/-", mi^p^} = 5xy5,,5ab5ap, {^r, = {^^^ = o, (5.3) 

where (a, b) and (a, /?), which have been suppressed in eq. ( |5.2| ), are color, flavor 
and Dirac indices, respectively. Of course, we have again replaced the classical link 
variables Ux,fj. by quantum link operators Ux,^- The generator of an SU{N) gauge 
transformation now takes the form 

Gx = E(^-A,M + Lx,^.) + ^lA^,, (5.4) 

and it is again straightforward to show that H commutes with Gx for all x. 

The same universality arguments that were used before now suggest that the 
effective action of the corresponding 5-d gauge theory (with = 0) takes the form 

Slip, ^, A^] = dx, J d^x {^[TrF^^F^, + ^Tid^A,d,A^] 

+ i^bM^^ + 9^.)+mq + ^^-i^^^]^|J}. (5.5) 

The "velocity of light" c', which characterizes the propagation of a quark in the fifth 
direction, is in general different from the corresponding quantity c for the gluons, 
because the quantum link formulation has no symmetry between the four physical 
space-time directions and the extra fifth direction. This is no problem, because we 
are only interested in the 4-d physics after dimensional reduction. 



17 



To ensure the proper dimensional reduction of the quarks, their boundary condi- 
tions in the fifth direction must be chosen appropriately. The standard antiperiodic 
boundary conditions, which are dictated by thermodynamics in the Euclidean time 
direction, would lead to Matsubara modes, ps = 271(77,5 + ^)//?, which would limit the 
physical correlation length of the dimensionally reduced fermion to 0{(3c'). The con- 
finement physics of the induced 4-d gluon theory, on the other hand, takes place at 
a correlation length which is growing exponentially with (3. In fact, (5 plays the role 
of the lattice spacing of the dimensionally reduced theory. Quarks with antiperiodic 
boundary conditions in the fifth direction would hence remain at the cut-off and the 
dimensionally reduced theory would still be a Yang-Mills theory without quarks. 
Once this problem is understood, one possible solution is obvious. One can simply 
choose periodic boundary conditions for the quarks in the fifth direction. This gives 
rise to a Matsubara mode, ps = 0, that survives dimensional reduction. Since the 
extent of the fifth direction has nothing to do with the inverse temperature (which 
is the extent of the Euclidean time direction), one could indeed choose the boundary 
condition in this way. 

Fortunately, we can do better. The above scenario with periodic boundary condi- 
tions for the quarks would suffer from the same fine-tuning problem as the original 
Wilson fermion method. The bare quark mass would have to be adjusted very 
carefully in order to reach the chiral limit. In practice this is a great problem in 
numerical simulations. This problem has been solved very elegantly in Shamir's 
variant of Kaplan's fermion proposal 0]. Kaplan studied the physics of a 5-d 
system of fermions, which is always vector-like, coupled to a 4-d domain wall that 
manifests itself as a topological defect. The key observation is that under these 
conditions a zero mode of the 5-d Dirac operator appears as a bound state localized 
on the domain wall. From the point of view of the 4-d domain wall, the zero mode 
represents a massless chiral fermion. The original idea was to construct lattice chiral 
gauge theories in this way. Shamir has pointed out that the same mechanism can 
solve the lattice fine-tuning problem of the bare fermion mass in vector-like theories 
including QCD. He also suggested a variant of Kaplan's method that has several 
technical advantages and that turns out to fit very naturally with the construction 
of quantum link QCD. In quantum link models we already have a fifth direction for 
reasons totally unrelated to the chiral symmetry of fermions. We will now follow 
Shamir's proposal and use the fifth direction to solve the fine-tuning problem that 
we would have with periodic boundary conditions for the quarks. 

The essential technical simplification compared to Kaplan's original proposal is 
that one now works with a 5-d slab of finite size (3 with open boundary conditions 
for the fermions at the two sides. This geometry limits one to vector-like theories, 
because now there are two zero modes — one at each boundary — which correspond 
to one left and one right-handed fermion in four dimensions. This set-up fits natu- 
rally with our construction of quantum link QCD. In particular, the evolution of the 
system in the fifth direction is still governed by the Hamilton operator of eq.( p.2D . 
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The only (but important) difference to Wilson's fermion method is that now r < 0. 
Of course, one could also obtain a left and a right-handed fermion by using a domain 
wall and an anti-wall with otherwise periodic boundary conditions. In that case the 
Hamiltonian of eg. (|5.2|) would have to be modified in an xs-dependent way. Shamir's 
method is more economical and concentrates on the essential topological aspects, 
which are encoded in the boundary conditions for the fermions in the fifth direc- 
tion. It is important that in Shamir's construction one also puts = 0. There are 
some minor differences between the implementations of the method in the standard 
formulation of lattice gauge theory and in quantum link models. In the standard 
formulation one works with a 4-d gauge field, which is constant in the fifth direction. 
In quantum link QCD this is not possible, because the non-trivial dynamics in the 
fifth direction turns the discrete states of quantum links into the continuous degrees 
of freedom of physical gluons. However, it is still true that the physical gluon field is 
essentially constant in the fifth direction, because its correlation length grows expo- 
nentially with (3. This is important for the generation of the fermionic zero modes 
at the two sides of the 5-d slab. 

The partition function of the theory with open boundary conditions for the 
quarks and with periodic boundary conditions for the gluons is simply given by 

Z = Tr(0|exp(-/3i/)|0). (5.6) 

Here the trace extends only over the gluonic Hilbert space of the quantum link model, 
thus implementing periodic boundary conditions for the gluons. We decompose the 
quark spinor into left and right-handed components 



= ( ;f 1 • (5.7) 



Open boundary conditions for the fermions are realized by taking the expectation 
value of exp{—(3H) in the Fock state |0), which is annihilated by all right-handed 



"^Rx and by all left-handed "^lx HlJ. As a result, there are no left-handed quarks 
at the boundary at Xs = and there are no right-handed quarks at the boundary 
at X5 = f3. Of course, unlike periodic or antiperiodic boundary conditions, open 
boundary conditions for the fermions break translation invariance in the fifth direc- 
tion. Through the interaction between quarks and gluons, this breaking also affects 
the gluonic sector. This is no problem, because we are only interested in the 4-d 
physics after dimensional reduction. For that it is essential that both quarks and 
gluons have zero modes, which is indeed the case with the boundary conditions from 
above. 



Of course, as we have argued before, the gluonic correlation length is not truly 
infinite as long as (3 is finite, but — due to confinement — it is exponentially large. 
As we will see now, the same is true for the quarks, but for a totally different reason. 
In fact, already free quarks pick up an exponentially small mass due to tunneling 
between the two boundaries, which mixes left and right-handed states, and thus 
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breaks chiral symmetry explicitly. To understand this, let us consider the actual 
Hamiltonian Ti-lp), which describes the evolution of free quarks with spatial lattice 
momentum p in time (rather than H, which describes the evolution of the system 
in the fifth Euclidean direction) 

nip) = 74[7555 + M+W{2 sin - 2 E 7. sinp,]. (5.8) 

i i 

The wave function of a stationary quark state with energy E{p) is determined by 

nimip,x,) = E{p}^{p,x,). (5.9) 

Let us consider states with lattice momenta = 0, vr. The state with p = 
describes a physical quark at rest, while the other states correspond to doubler 
fermions. Denoting the number of non-zero momentum components = vr by 

^ = Esin'|' (5-10) 

i 

the above Dirac equation ( ^.9|) takes the form 

-d^-^Lip, X5) + (M + 2nr)^L{p, X5) = E(p)^'r(p, Xs), 

d5^R{p, X5) + {M + 2nr)^^j(p, X5) = E{p}^l{p, x^). (5.11) 

Due to the boundary conditions introduced before, we must solve these equations 
with ^'l(p, 0) = ^'r(p, /3) = 0. Inserting the ansatz 

\E'l(p, X5) = Asinhaxs, '$r{p,X5) = ±Asmh.a{x5 - (3), (5.12) 

immediately gives 

a = (M + 2nr) tanha/3, M + 2nr = ±E{p) cosh a/3. (5.13) 

This equation has a normalizable solution only if M + 2nr > 0. In that case, for 
large (3 this implies a = ±(M + 2nr) and hence 

E{p) = ±2(M + 2nr) exp(-(M + 2nr)p). (5.14) 

The important observation is that the bare mass parameter M is not the physical 
quark mass, although we have considered an unrenormalized free theory. Let us first 
consider the doubler fermions, which are characterized by n > 0. It is then possible 
to choose r such that M + 2nr < 0, for which no normalizable solution exists. Thus, 
for r < — M/2, i.e. for a sufficiently strong Wilson-term with an unconventional 
sign, the doubler fermions are removed from the physical spectrum. The mass of 
the physical fermion (characterized by n = 0) is 

^ = \E{6)\ = 2Mexp{-Mp), (5.15) 
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which is exponentially small in (3. This situation is illustrated in fig. 2. The above 
result suggests how the fine-tuning problem of the fermion mass can be avoided. 
The confinement physics of quantum link QCD in the chiral limit takes place at a 
length scale 

1 , 24n'P , 

m°^^"P^(lliV-2iV,)e2)' ^''^'^ 

which is determined by the 1-loop coefficient of the /3-function of QCD with Nf 
massless quarks and by the 5-d gauge coupling e. As long as one chooses 

247r2 

the chiral limit is reached automatically when one approaches the continuum limit 
by making f3 large. For a given value of r one is limited by M < — 2r (note that 
r < 0). On the other hand, one can always choose J (and thus e^) such that the 
above inequality is satisfied. 

Of course, we also want to be able to work at non-zero quark masses. Following 
Shamir, we do this by modifying the boundary conditions for the quarks in the fifth 
direction. Instead of using iPl{x, x^, 0) = ipR^x, x^, /3) = we now demand 

ipiix, X4, 0) = -^^i-l^' ^4, f3), ^Pr{x, X4, 13) = ^4, 0), (5.18) 

where rriq is a mass parameter. This reduces to the previous boundary condition 
for rriq = 0, while it corresponds to antiperiodic boundary conditions for rriq = 2M, 
and to periodic boundary conditions for mg = —2M. Solving the above equations 
with the new boundary condition indeed yields physical quarks of mass in the 



continuum limit /? — oo as long as nig <^ M. It has been shown in ref.|2I| that in 
the interacting theory niq is only multiplicatively renormalized. On the level of the 
partition function 

Z = Ti[exp{-/3H)0{mg)] (5.19) 
the new boundary condition manifests itself by a mass- dependent operator 

0{mg) = YlC^nA + ^'^L'^R^){^lA + ^^L^L.), (5.20) 



which was constructed in ref.[^. Note that in eq.( |5.19| ) the trace is both over 
the gluonic and over the fermionic Hilbert space. In the chiral limit nig = the 
operator 0{mg) reduces to a projection operator |0)(0| on the Fock state introduced 
before. For rrig = 2M, i.e. for antiperiodic boundary conditions, the operator 0{mg) 
becomes the unit operator and the partition function reduces to the expression well- 
known from thermodynamics. 



Like any massive fermion action in odd dimensions, the effective action of eq. ( 
is not parity symmetric. Hence, there is the potential danger that the correspond- 
ing symmetry breaking terms make their way into the dimensionally reduced 4-d 
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theory. Indeed the topological 6'-vacuum term violates parity in four dimensions. 
For the moment we are interested in QCD with 6 = and therefore in a parity 
symmetric theory. Later, we will ask how to incorporate 6. We will now show that 
the dimensionally reduced theory is indeed parity invariant. For simplicity this dis- 
cussion will be based on the low energy continuum effective action of eq.( |5.5|) . It is, 
however, straightforward to apply the same arguments to the lattice action in the 
path integral that one obtains from the Hamiltonian of eg. ( ^.21) . The 5-d effective 
action is not invariant under parity as it is usually defined in five dimensions. Since 
we are interested in 4-d physics, this definition of parity is, however, not relevant 
for us. Instead, we now define another transformation that leaves the 5-d action 
invariant and reduces to ordinary parity in the dimensionally reduced 4-d theory. 
Let us consider 

ip'{x, X4, xs) = 74V'(-^, a;4, /? - X5), ij'{x, X4, X5) = ^(-f , X4, (3 - ^5)74, 
A[{x, X4, X5) = -Ai{-x, X4, /3 - xs), A'^{x, X4, X5) = A4(-f, X4, P -X5). 

(5.21) 

It is straightforward to show that this is indeed a symmetry of the 5-d effective 
action. Note that we have changed X5 to f] — X5. This change is invisible from the 
point of view of the dimensionally reduced 4-d theory. In fact, in the 4-d world 
the above transformation reduces to ordinary parity. However, exchanging x^ and 
(3 — X5 is essential as far as the 5-d theory is concerned. This operation exchanges 
the two domain walls and hence the left and right-handed fermions that are bound 
to them. This is indeed what parity is supposed to do. 

A non-perturbative formulation of QCD would be incomplete without a discus- 
sion of the vacuum angle 6. Although in the real world this parameter is indis- 
tinguishable from zero, it represents a possible parity breaking effect in the QCD 
Lagrangian. In the standard formulation of lattice QCD such a term could be added 
via the topological charge, for which a lattice regularized expression with adequate 



topological properties exists |2^. In quantum link QCD this would not be possible, 
because the discrete states of the quantum link operators can not directly encode 
topological properties. Here we propose to include exp(z6') as the determinant of the 



fermion mass matrix and hence via the boundary condition of eq.( ^.18|) . When one 



uses a fermion formulation which does not require fine-tuning in the chiral limit, this 
is the most natural thing to do, even in the standard formulation of lattice QCD. 



6 Glueballs, Mesons and Constituent Quarks 

In this section we express the Hamiltonian of the U{N) quantum link model with 
quarks in terms of color neutral operators. Gauge invariance requires that these 
operators are local bilinear combinations of rishons and quarks, which we refer to 
as glueballs, mesons and constituent quarks. All these objects — including the 
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constituent quarks — are bosons. In fact, the resulting representation of the theory 
can be viewed as a first step towards a bosonization of U (N) QCD. The determinant 
term in the SU{N) quantum link Hamiltonian would give rise to an additional color- 
singlet rishon-baryon consisting of rishons. For odd this object would hence 
be a fermion. To avoid complications related to these objects we limit ourselves to 
U (N) in this section. 

As in the Yang-Mills case, we represent the quantum link Hamiltonian of U {N) 
QCD in terms of rishons 



+ ^ E [(*Ic^,+m)757m(4+a -m*^+a) - (^l+A^c^+A -m)757/.(4,+,.*x)] 

X 

(6.1) 

Note that some of the terms in brackets now represent color neutral bosonic con- 
stituent quarks formed by a quark and a rishon located at the same end of a link. 
In fact, the Hamiltonian can be expressed as 

H —J ^x+jl, — fX, + V^X+jX-\-i>, — V,—lJ^X+V,+lJL, — U^X,-\-V,-\-IJ, 

X^IJL^V 

+ ^ E [<51,+m757m<3^+A -M - <5i+A -m757m<5x,+m] + E Tr(A475) 

^ x,it X 

+ ^ E [2Tr (M,75) - g;+^75g.+A,-M - Q1+a,-m75Q.,+/.] • (6.2) 



x,p. 



The traces are over flavor and Dirac indices. We have introduced 8^ = 64 glueball 
operators 

^x,±n,±u = E (^x,±A,±u^ (6-3) 

i 

which satisfy the local commutation relations of U (8) 

[^x,±iJ,,±iy, ^y,±p,±a] — 3xy{3±v,±p^x,±n,±a — 3±iJ.,±a^x,±p,±i^) ■ (6.4) 

We have also defined {ANfY meson operators 

i 

which carry flavor and Dirac indices and which generate the algebra of U {4:Nf) 

[Mr'^ M^'^'^'] = 5xy(5a/3jM^'''' - 5ad5asM'/'^^). (6.6) 
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While glueballs and mesons are not related via their commutation relations, i.e. 

[$x.,±m,±.,M;"^^] =0, (6.7) 
they are both related with the constituent quark operators 

Qr±M = E4!±,^r, (6.8) 

i 

via the commutation relations 

[*.,±M,±- Qr±p] = WT'^^ QV,±,] = S.ySaJa,QT,±,- (6.9) 

Finally, the commutation relations of the constituent quark operators take the form 

[QZ±,, Ql%u] = [QTl„ Ql%.] = 0. (6.10) 

Thus, the inclusion of the constituent quark operators completes the site-based 
algebra of U{8) ® U{ANf) to ?7(8 + 4A^/). We expect this representation of the 
theory to be useful in the investigation of the large N limit and in attempts to 
bosonize quantum link QCD. 



7 Conclusions 



We have constructed a new lattice formulation of QCD in the framework of quantum 
link models. In these theories the link-based Hilbert space of the gluons is finite. 
This is possible because the theory is constructed with a fifth Euclidean direction 
that ultimately disappears via dimensional reduction, although in the continuum 
limit its extent (3 diverges in units of the lattice spacing of the quantum link model. 
On the other hand, 13 itself plays the role of the lattice spacing of another effective 4- 
d lattice theory, whose non-perturbatively generated length scales are exponentially 
large in [3 (see fig. 3). 

An important question is whether SU{N) quantum link models on a 4-d space- 
time lattice possess a continuum limit in the universality class of the Coulomb phase 
of 5-d non-Abelian gauge theories. This is essential for their proper dimensional 
reduction to ordinary 4-d Yang-Mills theories. One way to investigate this is to 
study the classical equations of motion of quantum link models. As will be shown in 
a forthcoming publication [^], they indeed resemble those of 5-d Yang-Mills theories 
with = 0. This shows that in the classical limit — i.e. when one works with 
a large representation of SU (2iV) — dimensional reduction does indeed take place. 
Whether dimensional reduction also occurs with the 20-dimensional representation 
of SU (6) and hence whether the SU{3) quantum link model introduced here provides 
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a viable cut-off scheme for tlie continuum field theory of QCD, is an issue that can 
only be decided in numerical simulations. This situation is analogous to the one 
in quantum spin models. Based on the classical equations of motion one finds 
that for large spin the 2-d antiferromagnetic quantum Heisenberg model has an 
ordered ground state and that the corresponding Goldstone bosons — in that case 
antiferromagnetic magnons or spin-waves — have a relativistic dispersion relation. 
Still, it required high-precision numerical simulations [ll| to be sure that a 



staggered magnetization is spontaneously generated even in the extreme quantum 
limit of spin 1/2. Only this finally justifies the use of a 3-d relativistic low-energy 
chiral Lagrangian with 0(3) symmetry to describe the dynamics of the spin 1/2 



quantum Heisenberg model at low energies ||2^. Similarly, at present we cannot 
guarantee that the low-energy excitations of a 4-d SU (3) invariant quantum link 
model formulated with the 20- dimensional representation of SU (6) are 5-d massless 
gluons. However, this should at least be the case in the classical limit, i.e. when we 
work with a large representation of SU{6). 

It is remarkable that the introduction of a fifth Euclidean direction of finite 
extent also naturally protects the chiral symmetries of lattice fermions, when open 
boundary conditions are used for the quarks at the two 4-d boundaries of the 5-d 
slab (see fig. 2). This set-up — first proposed by Shamir as a variant of Kaplan's 
fermion method — fits very naturally with quantum link models. Due to tunneling 
between the two boundaries, a fermionic correlation length is generated, which is 
also exponentially large in /3. Quantum link QCD treats quarks and gluons on an 
equal footing. They both have finite Hilbert spaces (per unit volume), both require 
a fifth Euclidean direction that finally disappears via dimensional reduction and 
they both have correlation lengths that are exponential in (3 without requiring fine- 
tuning. We believe that these properties point to an interesting connection between 
fermions and gauge fields in quantum link models. This may help in attempting to 
formulate chiral or supersymmetric gauge theories on the lattice. 

Due to the fact that their gluonic Hilbert space is finite, quantum link models 
can be expressed in terms of rishons, which are fermionic constituents of the gluons. 
The rishon dynamics can probably not be formulated as a continuum field theory. 
Therefore we consider the rishons as lattice artifacts in our way to formulate QCD. 
The gluons in quantum link models are composites of rishons, just like antiferro- 
magnetic magnons can be viewed as composites of electrons that are hopping on a 
lattice. Of course, in that case we know that the electrons are real particles, which 
are hopping on a real crystal lattice. If our world is a thin five- dimensional slab with 
a lattice structure at extremely small distances, the rishons in quantum link models 
could be real particles. We do not want to speculate about that possibility and thus 
we consider them as mathematical objects only. However, the rishons may turn out 
to be of great importance for solving QCD in the large limit. We consider this 
possibility the most promising analytic aspect of the new formulation. 

At finite we do not see a way to solve quantum link models analytically. How- 
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ever, from a computational point of view quantum link QCD is still very attractive, 
in particular because the new framework may allow one to construct very efficient 
cluster algorithms. The five-dimensional set-up is not as inconvenient as it may 
seem, because — due to the discrete nature of the Hilbert space — the 5-d path 
integral for quantum link models can be formulated and — most important — also 
be simulated directly in the continuum of the fifth Euclidean dimension. In fact, 
even in the standard formulation of lattice QCD it may turn out to be advantageous 
to work in five dimensions in order to control the chiral properties of the quarks. 

It remains to be seen if quantum link models provide a more efficient formulation 
of QCD than standard lattice gauge theory. In any case, quantum link models allow 
us to attack the long-standing QCD problem from a different perspective. 
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